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The dynamic behavior of the entanglement for two two-level atoms coupled to a common lossy cavity is 
studied. We find that the speed of disentanglement is a decreasing (increasing) function of the damping rate 
of the cavity for on/near (far-off) resonant couplings. The quantitative explanations for these phenomena are 
given, and further, it is shown that they are related to the quantum Zeno and anti-Zeno effect induced by the 
non-Markovian environment. 
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I. INTRODUCTION 

Quantum entanglement plays an essential role in quantum 
information science llfl. However, entanglement is very frag- 
ile due to the influence of the environment, which is recog- 
nized as the major obstacle for the final application of quan- 
tum information processing. For example, a phenomenon 
called "entanglement sudden death" 10. S H] shows that two 
initially entangled qubits may disentangle completely in a fi- 
nite time due to spontaneous emission. So in order to keep the 
atomic entanglement for a long time, spontaneous emission 
should be suppressed. Several ways were proposed for this 
purpose. One way widely applied is to place the qubits in a 
structured environment, say, microcavity J5l |6[] or in the pho- 
tonic band gap of photonic crystals 01, such that the qubits 
are separated from the environment. Another way is to dy- 
namically control the coupling between the system and the 
environment, by, e.g., quantum Zeno effect ML In prac tice, 
people may combine the two ways ||^[^[l^[lli[H[l3i[l4ll. 
In fact, the two methods have something in common. 

In this paper, we investigate the entanglement dynamics of 
two two-level atoms in a common cavity, a non-Markovian 
environment. We find that when the atoms are on or near res- 
onant with the cavity, the speed of the disentanglement de- 
creases as the quality factor of the cavity decreases, and when 
the atoms are far-off-resonant with the cavity, the speed of the 
disentanglement increases as the quality factor of the cavity 
decreases. These phenomena can be related to q uantum Zeno 
OSBGl and anti-Zeno effect 03 El EH. 



II. THEORETICAL MODEL 

Consider two spontaneously emitting two-level atoms with 
a common zero-temperature bosonic reservoir, and in the 



Lamb-Dicke limit lfl7ll . the dipole-dipole interaction is neg- 
ligible. Under the rotating wave approximation (RWA), the 
Hamiltonian of this composite system plus the reservoir is 
given by (ti = 1) flEOJ] : 

H = H +H int , (1) 

with 

Hq — a>icr+cr_ + afycr+cr^ + I dcj^kb' \ojk)b((Dk), 



H int = (fficr^ + ff 2 o~ ( + 2) ) I dto k g(a) k )b(aj k ) + h.c. 

xj — OO 

Here, <r± and cjj are the inversion operators and transition 
frequency of the jth qubit (j = 1,2), and b(o)k), b^(a>k) are the 
annihilation and creation operators of the field mode of the 
reservoir. The mode index k contains several variables which 
are two orthogonal polarization indices and the propagation 
vector k. To measure the coupling strength of the atoms to 
the cavity mode determined by the atom's relative position 
in the cavity, we introduce the dimensionless constant aj JH]. 
In the following discussion, we introduce vacuum Rabi fre- 
quency R = W(a 2 x + q-j) 1 ^ 2 and relative coupling strengthes 

rj = aj(al+alr V2 U = h2). 
For an initial state of the form 

|0<O)> = (ci O |e>ilg>2+c 2 ol£>ile> 2 )|0> £ , 

dto k b [ (LOk)b(a> k ) + 

cr+V^ + cr+V?\ the time evolution of the total sys- 
tem is confined to the subspace spanned by the bases 
{\e)i\g) 2 \0)E,\g)i\e)2\0)EAg)i\8)2\W E }- 

\m) = c^e-™" \e\ \g) 2 \0} E + c 2 (t)e-^' \g\ \e) 2 \0) E 

/-*oo 

+ dcokC^We-^'lgWg^lWe, (2) 
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where is the state of the reservoir with only one exciton 
in the kth mode. Here, we consider the case in which the two 
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atoms have the same Bohr frequency, i.e., a>i = io 2 — (jl>o- 
Substituting Eq. (f2]i into Schrodinger's equation and eliminat- 
ing the coefficients c Mt (f), one has 



61(f) 



f 

Jo 



dhf(t - ti)aii[aici(ti) + a 2 c 2 (h)], 



(3a) 
(3b) 



6 2 (t) = - I dhf(t - ti)a 2 [aici(ti) + a 2 c 2 (ti)], 
Jo 

where the correlation function takes the form: 

da> k J(a) k )e- i ^-'*>*-V. 

As discussed in [8], there is a subradiant state, |^_) = 
r 2 \e)i \g) 2 - r\ \g) 1 \e) 2 , which does not decay in time, and the 
only relevant time evolution is the superradiant state \tff+) - 
r\ \e)i \g) 2 + r 2 \g) 1 \e) 2 . 

For the case of two atoms interacting with a cavity field in 
presence of cavity losses, the spectral density function takes 
the form 



J(co k ) = W 2 A/n[(co k - u c ) 2 + A 2 }, 



(4) 



where W is the transition strength, (±> c is the center of the spec- 
trum, and 2A is the full width at half maximum (FWHM) 
of the spectral function. By employing Fourier transform 
and residue theorem, we get the explicit form f(t — ti) = 
yf2 e -A\t-h\^ wnere tjjg quantity l/Ais the reservoir correlation 
time. 

Using Laplace transform, we get the solutions of Eqs. < f3ab 
and Ob}: 



ci(f) = r 2 (3- 
c 2 (f) = -nP- 



nP+e(f), 
r 2 p + s(t), 



(5a) 
(5b) 



where p ± = {ip ± \ \{r(f)),s(f) = (s++A + i6)e s+ '/(s + -s-)-(s- + 
A + iS)e s -' l(s+ - s~), and s± are the roots of the equation for 
s: s 2 + (A + i6)s + R 2 — 0, where 6 — u> c - u>o is the detuning. 
In the {|e>i \e) 2 , \e) l \g) 2 , \g) x \e) 2 , \g) t \g) 2 ] basis, the reduced 
density matrix of the two atoms is given by: 



Pa(f) 



( 

| Cl (f)| 2 

c 2 (t)c\(f) 






ci(t)c* 2 (t) 
\ci(t)\ 2 








1 - \ci(f)\ 2 - \c 2 (t)\ 2 



(6) 



The entanglement of the two atoms can be evaluated by con- 
currence C(t) IU9I1 - For p a [Eq. ©], its concurrence can be 
derived from 11911. as 



C(0 = 2|ci(/)c;(0|=2|ci(0l|ca(0l. 



III. NUMERICAL RESULTS 



(7) 



We focus on the concurrence as a function of time t in 
weak-coupling regime, W < A/2. In this regime, the concur- 
rence of the two atoms undergoes nearly irreversible exponen- 
tial decay. Similar behaviors mentioned below take place for 
strong-coupling regime in the time scale of Rabi oscillation. 




(a) 




(b) 



FIG. 1: (Color online) Time evolution of the concurrence, with (a) 
<5 = and (b) S = 20, both with the initial state being the maximal 
entanglement state \(p+), W = 1, and r\ = r%, for the cases of (i) 
A = 5 (solid blue curve), (ii) A = 8 (dotted green curve), (iii) A = 10 
(dashed red curve) . 



We compare the entanglement dynamics of the two 
atoms initially in the maximal entanglement states \tp + ) = 
"^d e )t 1^)2 + l e )2) f° r three different values of full width 
at half maximum (FWHM) of the spectral function, namely, 
2^=10,16,20. 

As mentioned in 12011 . it is hard to exactly control the posi- 
tion of the atom in the optical cavity. But through numerical 
stimulation, we find that different values of r\ show qualita- 
tively similar behaviors. So, for simplicity, we focus on the 
case of equal coupling parameters, i.e., r\ — r 2 — V2/2. 

As in FIG. |TJa), the atoms are on resonance with the center 
of the spectrum, (5 = 0. The concurrence decreases mono- 
tonic ally down to zero at the beginning. An interesting phe- 
nomenon is that the speed of disentanglement decreases as A 
increases. Similar behavior happens when the atoms are near 
resonance with the center of the spectrum, 6 <K A. In fact, this 
phenomenon is related to the environment induced quantum 
Zeno effect iTToL Il5[ [l6ll . However, When the atoms are far 
off-resonant with the center of the spectrum, 8 » A, for ex- 
ample, as shown in FIG. [Tib), where we choose 6 = 20, the 
speed of disentanglement increases as A increases. This latter 
phenomenon is related to the anti-Zeno effect 111 l[ll2lll3lll411 . 

In fact, as mentioned in fall , if the coupling strengths of 
the two atoms to the field are different and the dipole-dipole 
interaction is not negligible, there are no asymptotic entangle- 
ment, which means that even \ifrJ) will disentangle completely. 
In these cases, our numerical stimulation shows that similar 
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FIG. 2: (Color online) The spectrum of the field in the cavity, /(A) = 
W 2 A/n(A 2 + A 2 ), W = 1, for the cases of (i) A = 5 (solid blue curve), 
(ii) A = 8 (dotted green curve), (iii) A = 10 (dashed red curve) . 



phenomena mentioned above happen. 

We can give an intuitive explanation for these phenomena. 
As we can see in FIG. |H the center part of the spectrum de- 
creases monotonically as A increases, while the parts which 
are far from the center increases as A increases. We can 
prove that the short-time behavior of the disentanglement is 
determined by the modes of the spectrum which are on res- 
onance with the atoms: the speed of the disentanglement de- 
creases (increases) as the density of these modes decreases 
(increases). 

Similar to (MCI, we define P(t) = |e(f)| 2 = e~ w , where 
R is the effective decay rate. From Eqs. d5ab . ( Bbl and (0, we 
can see that P(t) and R can describe the concurrence to some 
extent. For short-time behavior, in the first-order approxima- 
tion, one yields 1 1311 : 



dcjJ(cj)F(aj), 

OO 

? . o (CD - COq\ 



(8a) 
(8b) 



Since F(a>) is a sharply steep function of a> around coq, if the 
width of F(u>) is much smaller than that of J{oS), R is mainly 
determined by J{(jJq). In this paper, the form of J(a>) is given 
by Eq. (0J. Then a detailed analysis shows that if A is larger 
than cl>o - a> c , J(a>o) is a monotonic decreasing function with 
respect to A, otherwise, it is a monotonic increasing function 
with respect to A. Since R oc J(cjq), these results also go to R. 



IV. ENVIRONMENT INDUCED QUANTUM ZENO 
EFFECT AND ANTI-ZENO EFFECT 



In order to relate these phenomena to environment induced 
quantum Zeno and anti-Zeno effect, we make a transform on 
the original Hamiltonian [Eq. (Q3] as in ll2~ll l22tl : 



a = (A/n) 



1/2 



f 



d(x>k- 



b(cj k ) 



iA' 



(9) 



c(A) = (A/n)P.V.{ daj, 



b(a> k ) 



' (a> k - A)(tu k -co c - iA) J 



we get 



(x>k — We 

du> k —b(u>k), 

OO 10k - lO C — lA 



H — Z/q + Hj n f ~t~ H ' damping i 



J4' - , , rr^lr^ 



(D_(l) 



, , , ^(2)^(2) 
+ Cl>2(T + <X_ 



H' mt = W[(aKT^ + a 2 oi 2) )a + (a { cr m + a 2 (r {2) )a\ 



(10) 



(ID 



(12) 



(13) 



Hdamping = I <iAAc + (A)c(A) 

U — OO 

,iA[fl t c(A) + flc 1 "(A)], (14) 

OO 

where a 1 and a are the creation and annihilation operators for 
the discrete quasimode, with the frequency co c , and c^(A), c(A) 
are those for the continuum quasimodes with the frequency A. 

It can be seen, from the above Hamiltonian, that the com- 
posite two-atom system only interacts with one discrete mode 
and their coupling coefficient is just the transition strength 
W. The discrete mode interacts with the continuum modes 
and their coupling strength only contains the width (A) of the 
Lorentzian spectral density which is a constant. This means 
that if we let the two atoms and the discrete mode be a new 
system and the continuum modes be the reservoir, the behav- 
ior of the new system is exactly Markovian. As discussed in 
1I22I1 . the discrete mode is appropriate for describing the elec- 
tromagnetic field inside the cavity, and the continuum modes 
are appropriate for describing the field outside the cavity. 

In terms of quantum measurement theory |23[ l24l l25ll . the 
continuum quasimodes c(A) represent the environment, the 
discrete quasimode a represent the meter, and the two atoms 
are the system. By interpreting the environment as a photo- 
electron counter, the coupling strength of interaction between 
the discrete quasimode and the continuum quasimodes relates 
to the rate of photon counting A. That means the environment 
acts as an observer of the pointer states of discrete quasimode. 
Although, in our model, |l) fl and |0) fl are not perfect pointer 
states for discrete quasimode, the observation by the environ- 
ment destroys the coherence between them. We can interpret 
t — A~ as the effective measurement interval by the environ- 
ment. Assume that the two atoms and the discrete quasimode 
are in the state |(/?+) 12 |0) fl in the beginning. For simplicity, 
we choose a\ = then because of the interaction between 
the atoms and the quasimode, at a subsequent time r, they 
evolve to the state |<^+)i2 |0) fl + WT\gg) 12 \l) a , and the prob- 
ability of the transition \ifi+) l2 \0) a — > \gg)i 2 |l) fl is (Wt) 2 . In 
the meantime, the coherence of this superposition is destroyed 
by the environment. Then the probability of the transition 
after a time At » t is P * (Wr) 2 Af/r = W 2 tAl So, we 
can see that this probability is proportion to A -1 , It explains 
the phenomenon that the speed of disentanglement decreases 
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as the strength of damping increases. In terms of quantum 
measurement theory, this means that the observation made by 
the environment suppress the disentanglement, which is just 
quantum Zeno effect. When ujq - a> c is much larger than A, 
that means the frequency of the observation is very small, and 
in the language in lfl2[ [A . we can say that the observations 
are made in the anti-Zeno regime, so the observations acceler- 
ate the disentanglement. A similar discussion about the Zeno 
and anti-Zeno effect for nonresonant systems is given in lfl4ll . 
where the observations are not made by the environment. 

In conclusion, we extend the study of the entanglement dy- 
namics of two atoms in a common cavity. We find that the 
speed of the disentanglement of the two atoms is a decreasing 



(increasing) function of the damping rate of the cavity when 
the atoms are on/near resonance (far off resonance) with the 
center of the cavity modes. We give a quantitative explana- 
tion for these phenomena, and relate them to quantum Zeno 
and anti-Zeno effect induced by the environment. These re- 
sults are helpful for understanding the related experimental 
phenomena and for the practical engineering of entanglement 
in the future. 
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